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Abstract 

We present a reciprocal space technique for the calculation of the Coulomb integral in two dimensions 
in systems with reduced periodicity, i.e., finite systems, or systems that are periodic only in one dimension. 
The technique consists in cutting off the long-range part of the interaction by modifying the expression 
for the Coulomb operator in reciprocal space. The physical result amounts in an effective screening of 
the spurious interactions originated by the presence of ghost periodic replicas of the system. This work 
extends a previous report [C. A. Rozzi et al, Phys. Rev. B 73, 205119 (2006)], where three-dimensional 
systems were considered. We show that the use of the cutoffs dramatically enhances the accuracy of the 
calculations, and it allows to describe two-dimensional systems of reduced periodicity with substantially 
less computational effort. In particular, we consider quantum-dot arrays having potential applications in 
quantum information technology. 

PACS numbers: 71.15.Mb, 73.21.La, 71.10.Ca 
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I. INTRODUCTION 



The interest in low-dimensional electronic structures has increased steadily during the past 
few decades. This is mostly due to breakthroughs in semiconductor technology in the 1970s 
and early 1980s. At present, low-dimensional systems form a significant portion of the whole 
field of condensed matter physics. Some examples are layered semiconductor devices such as 
metal-oxide- semiconductor field effect transistors, quantum Hall systems, spintronic devices, and 
quantum wells, wires, and dots 1 (QDs). 

Following the advances in constructing techniques for single two-dimensional (2D) QDs with 
tunable atom-like properties, it has become possible to couple QDs to form artificial "molecules". 
Coupled QDs have significant potential for solid-state quantum computation through, e.g., coher- 
ent manipulation of spins. 2 Furthermore, extended lattices or arrays of QDs have been fabricated, 3 
which, in addition to the proposed applications in quantum information, 4 show interesting mag- 
netic phase transitions 5 ' 6 which may be exploited in quantum transport and spintronics. 

From the theoretical point of view, dealing with systems of arbitrary periodic dimensionality 
may be complicated. In the simplest case of a fully periodic lattice of elements, periodic boundary 
conditions are applied at every cell border, and Bloch's theorem describes the discrete-translation 
invariant form of the orbitals. On the other hand, in all the cases in which the system has reduced 
periodicity, the use of a supercell with periodic boundary conditions in all the directions becomes 
problematic. In fact, the response function of a periodic lattice is generally very different from the 
response of a system with reduced periodicity (such as an isolated system, a chain, a slab, etc.) 
and the convergence of the fully periodic quantity to the reduced-periodic ones as a function of 
the supercell size is very slow. A large supercell is a numerical disadvantage, but it is specially 
necessary to avoid the influence of the periodic images if long range operators are used. 

The main issue here is indeed the computation of a long-range operator, i.e., the Hartree (or 
"Coulomb") potential: 



which is ubiquitous in Science, and which we study here in the context of 2D electronic structure 
calculations. In particular, we will exemplify our approach by utilizing density-functional theory 
(DFT), although the method we propose can also be useful in different fields. In Fourier space, the 
convolution integral (OQ) is transformed into a trivial product. This fact adds to the other undoubted 
advantages of the supercell approach, such as the natural inclusion of the periodic boundary con- 
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ditions, and the existence of very efficient fast Fourier transform algorithms. 

The attempt to retain these advantages, i.e., to compute the Hartree integral in reciprocal space, 
has led to the creation of several cutoff schemes for finite systems, whose main intent is to provide 
an effective truncated Coulomb interaction such that the system becomes unaware of the existence 
of its periodic replicas.-^ More recently, an exact scheme was proposed to achieve a broader 
goal, 10 namely to truncate the Coulomb interaction in a 3D system in the dimensions along which 
the system is confined, leaving it long-ranged in the dimensions in which the system is periodic. 

In this paper we show that a similar scheme can be drawn in a 2D space, allowing us to correct 
the spurious supercell effect when treating finite systems and one-dimensional (ID) chains. In 
particular we focus on the case of a single infinite chain of 2D few-electron QDs, that, in a classical 
supercell approach, would mistakenly appear as periodic in both directions, while it should be 
treated as a truly periodic system only along the x direction. We study the magnetic ground state 
of this system, and show that the use of the cutoff allows to speedup the calculations of ground 
state fundamental quantities such as the Fermi level. The technique is exact if the computation 
parameters are chosen judiciously. The possibility of isolating the chain from the replicas permits 
to get an insight about the role of inter-chain interaction in determining the conductive of insulating 
character of the chain. 

II. METHOD 

We follow closely the procedure described in Ref. Iioll . The Hartree integral is the convolution 
of the charge density n with the Coulomb interaction potential v. In 2D, 



We consider the charge density to be in a unit cell. This unit cell may then be replicated in 
both directions to fill all 2D space (the two dimensional periodic case, or "2D/2D"), replicated 
in only one direction ("2D/1D"), or not replicated at all (finite case, or "2D/0D"). If we move 
to reciprocal space, the unit cell is always replicated periodically in both directions, which is 
undesired in both the 2D/1D and 2D/0D cases. The Hartree integral in reciprocal space reduces to 
the simple multiplication: 




space 



(2) 



V(G)=h(G)v(G), 
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where G = (G x ,G y ) are the reciprocal vectors. The Fourier transform of v(x,y) can be readily 
computed: 

/.+00 /.+00 e i(G x x+G y y) 2% 

V {G x ,G y ) = jjxjjy^= = -. (4) 

This expression implies full 2D periodicity, and therefore it contains spurious terms if the period- 
icity is reduced. Our aim is to modify the expression of the Coulomb interaction v(G) — > v(G) to 
an effectively truncated interaction: 

{i if re© 
(5) 
ifr£z>, 

for some suitable region £> , such that it avoids the interaction of the "real" cells with the spurious 
replicas, while maintaining all interactions between points in the real cells. In order to achieve this 
goal, it will be necessary to increase the size of the original unit cell. 

For finite systems (2D/0D) the cutoff region can be conveniently defined as r < R, for some 
cutoff radius R. We can then easily perform the Fourier integral in polar coordinates: 

"2tz 



v 0D (G) = dr <*t>exp(iGrcos<|>) 
Jo Jo 

2n£drJ Q (Gr) = (2nR) X F 1 
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where 7o is the Bessel function of the first kind, and 1F2 is the generalized hypergeometric function. 
Note that the G = case is finite and continuous, and poses no difficulty: v 0D (G = 0) = 2%R. The 
value of R must be sufficient to contain all possible interactions in the original unit cell where the 
charge is contained; if we imagine it to be a square of radius L, then R = s/lL. If we now enlarge 
this cell (padding the density with zeros) to Ll = ( 1 + v / 2)L, the spurious replicas will not interact 
thanks to the interaction cutoff, and the Hartree integration will be exact. 

The 2D/1D case is more subtle. We assume the charge density to be contained in a strip defined 
by \y\ < R/2: the system is a chain of unit cells along the x axis. We define the cutoff region <D 
as \y\ < R - and the unit cell is also enlarged in the y direction to \y\ < R. It is easy to see that in 
this manner, the "ghost" replicas do not interact with the original chain, but the real interactions 
are preserved. The reciprocal space expression for the truncated Coulomb potential is: 

r+°° r+R Q i{G x x+G y y) 

v 1D (G x ,G > ,)= / dx dy (7) 

J -00 J-R ^x 2 +y 2 

= 4 [ R dycos(G y y)Ko(\G x y\), (8) 
Jo 
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where Kq is the modified Bessel function of the second kind. However, in this case, since 
lim Gj .^ + Ko(G x y) = +°°, the integral is undefined on the whole line G x = 0. 

In the fully periodic (2D/2D) case, we also have a singular point at G = [see Eq. ©]. If 
we assume charge neutrality, this singularity poses no problem. In Eq. © the term v(G = 0) 
multiplies n(G = 0), which is zero because this is in fact the charge neutrality condition. In the 
2D/1D case, we must estimate how these divergent terms affect Eq. ©, or, more precisely, the 
back Fourier transform 

y(r) = f d 2 Ge- iG - r n(G)v(G) . (9) 

In order to see how the infinities appear, we consider the integral in Eq. © for G x = 0, but using 
first a finite integration domain also in the x direction, — h < x < h. This integral is convergent. 
If we perform the integration and retain only the terms that do not vanish in the limit h — ► °° we 
obtain 

v 1D (0, G y ) « 41og(2/*) sm(G y R)/G y 

-4 f dycos(G y y)log(y). (10) 
Jo 

The first term, which we call v°°(0, G y ), diverges as h — > °°. However, it can easily be seen that 
it can be ignored if we assume charge neutrality. We perform the G y integration in Eq. © for 
G x = 0, considering only the v°°(0, G y ) term 

f dG y v°°{0,G y )n(0,G y )e- iG yy = 
41og(2A) f dx' [ m dy'n(x\y>) [ dG,^^^^ . (1 1) 

J J-R/2 J ' Gy 



Now we have: 



/ 



su^ e ^ y) = ^/2 « W-y\<R (12) 

Gy 



In the integral \y'\ < R/2, since the charge is contained in that region. This is the region of in- 
terest, and therefore we are also interested in looking at the potential only for \y\ < R/2. As a 
consequence, |/ — y\ < R, and we can conclude 



J dG v v oo (0,G > ,)n(0,G > ,)^ !GvV = 
r r R / 2 

2n\og(2h) dx' dy'n(x',y'). (13) 

J J-R/2 
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This is obviously zero if we assume charge neutrality. Therefore we have prooved that we can 
safely ignore the diverging terms, and retain only the regular ones. With this in mind, the case of 
G = is now trivial: v 1D (0,0) = -4R(logR - 1). 



HI. RESULTS 




To test the cutoff method, we consider ID QD arrays similar to those in Ref. Q6Q . Each rect- 
angular unit cell contains two QDs, and each QD has N electrons bound by a Gaussian positive 
background charge density. The total Gaussian background charge density has the form 



where r = (x, y), R = (na x ,0) with n = 0, 1,2, .. ., and r s = 2 is the average density at the center of 
the QD. Note the use of eff. a.u. throughout.— 

We solve the Kohn-Sham equations within spin-DFT on a 2D grid with and without the cutoff 
method described above. In the former case, the system is periodic only in x direction, whereas 
in the latter case it is periodic in both x and y directions. For the exchange and correlation we 
use the 2D local-spin density approximation (LSDA) with the parametrization of the correlation 
by Attaccalite et al. n This parametrization has been shown to be more accurate than the form of 
Tanatar and Ceperley 13 in the partially spin-polarized regime. 14 All the numerical calculations are 
done using the octopus^- code. 

In Fig. Q] we show the Fermi energy Ep of a QD chain with three electrons per dot as a function 
of the lattice constant a y perpendicular to the chain. The lattice constant in the x direction is fixed 
to a x = 5.06. Using the cutoff method leads to a very fast convergence, so that at a y = 3a x the 
Fermi energy is converged to six digits. In contrast, without cutoff, i.e., the system being periodic 
in both directions, a much larger supercell is needed in order to achieve comparable accuracy. This 
is due to the long-range Coulomb interaction between parallel QD chains. Therefore, for accurate 
calculations the cutoff scheme is essential in reducing the cell size and thus the computational cost. 

Next we focus our attention to the physical effects caused by multiple parallel QD chains in 
comparison with a single QD chain. In the former case no cutoff is used so that for a y = a x we 

have parallel chains located at y = na x , n = ±1, ±2, In the latter case (single chain) we use 

the cutoff method with a y = 3a x to guarantee a high precision (see Fig. [[]). 

In Fig. [2] we plot the band structure for a QD chain with three electrons per dot at a x = 5.06 
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Figure 1 : (Color online) Fermi energy of a quantum-dot chain with three electrons per dot as a function of 
the supercell size perpendicular to the chain. 

(left panel), corresponding to the system considered in Fig.Q] We find that the presence of parallel 
chains leads to a rather minor effect on the bands. The qualitative shape of the structure is very 
similar, and the shift in the bands due to the interchain effects is of the same order of magnitude 
as the shift in Ep. 

Fig. [2] can be also directly compared with the results of Ref. |]6Q. Overall, we find an excellent 
agreement, also in the case of a larger lattice constant a x = 8.55 (right panel). Here the spin 
degeneracy is lifted due to the exchange effect which leads to a magnetic ground state. 6 

Next we consider the interchain effects on phase transitions in QD chains. As noted by 
Karkkainen et al.,— the chains have a rich phase diagram with respect to the electron number 
N and the lattice constant a x . In general, increasing a x leads to insulating behavior and/or mag- 
netism, i.e., spin polarization, depending on N. In Fig. [3] we show the band structure at a x = 5 
for two electrons per dot, calculated with the cutoff procedure for a single chain and without the 
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Figure 2: (Color online) Band structure at lattice constants a x = 5.06 (left) and a x = 8.55 (right) for a 
quantum-dot chain with three electrons in each dot. The dotted lines in the left panel correspond to the 
result with multiple parallel chains at y = na x = na y , n = ±1, ±2, .... The straight lines correspond to the 
Fermi energies. 

cutoff for multiple parallel chains located at y = na x = na y , n = ±1, ±2, The single QD chain 

is clearly a metal with bands crossing the Fermi level. In contrast the presence of periodic replicas 
of the chain opens up a gap across the Fermi level so that the system becomes an insulator. The 
physical origin of the effect is the increased localization of single QDs due to Coulomb repulsion 
between the replicas of the chain. Similarly, the single QD chain (left panel in Fig. [3]) becomes an 
insulator at around a x = 6, when the dots are more isolated. 

Finally, we have also examined the interchain effects on the magnetism of QD wires. Generally, 
the spin polarization is relatively stable: For N = 1 ... 5 we find similar spin-polarization as a 
function of the lattice constant regardless of the presence of parallel QD chains. In all cases the 
magnetization is in a qualitative agreement with the results in Ref. [|6], minor differences arising 
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N = 2, a x = 5 




Figure 3: (Color online) Band structure at a lattice constant a x = 5 for quantum-dot chains with two electrons 
in each dot. The left and right panels show the results for a single and multiple chains, leading to to 
metallic and insulating band structures, respectively. In the latter case, the parallel chains are located at 

y = na x = na y , n = ± 1 , ±2, 

from the use of the LSDA parametrization by Attaccalite et al. instead of the one by Tanatar and 
Ceperley. 13 However, further studies are required regarding the validity of the LSDA in the high- 
correlation regime and/or in the presence of magnetic fields. These would also be ideal systems to 
test the 2D density functionals recently developed . 1 17 

IV. SUMMARY 

The reciprocal space is the natural venue for treating periodic systems. If the periodicity is not 
complete (the system is not periodic in all the space dimensions) or even absent, the computations 
require the use of a large supercell in order to avoid the spurious interactions due to "ghost" system 
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replicas; this is specially manifested in the computation of the Coulomb, or Hartree, integral, 
which is a trivial multiplication if we use plane waves. For 3D systems, it has been shown that the 
use of a screened Coulomb interaction greatly improves the accuracy in the calculation of ground 
states quantities, and substantially simplifies the evaluation of excited-state properties of reduced- 
periodicity systems. In this work we have shown that the same ideas can be applied to the case 
of the 2D electron gas; we have provided the relevant formulae for finite systems in 2D, and for 
systems that are periodic in only one dimension. 

Moreover, the corrective cutoffs are exact and rather straightforward to apply. We have demon- 
strated their effectiveness by computing band structures and Fermi energies of one-dimensional 
periodic arrays of quantum dots formed by the confinement of two-dimensional electron gas. We 
expect that our technique will be of great interest for studying these novel systems of "artificial 
molecules" and "crystals", although the ubiquity of the problem we have faced will probably make 
our technique relevant in other fields, too. 
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